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Abstract 

Using some simple toy models, we explore the nature of the brane-bulk interaction for cosmo- 
logical models with a large extra dimension. We are in particular interested in understanding the 
role of the bulk gravitons, which from the point of view of an observer on the brane will appear 
to generate dissipation and nonlocality, effects which cannot be incorporated into an effective 
(3+l)-dimensional Lagrangian field theoretic description. We explicitly work out the dynamics 
of several discrete systems consisting of a finite number of degrees of freedom on the boundary 
coupled to a (l+l)-dimensional field theory subject to a variety of wave equations. Systems both 
with and without time translation invariance are considered and moving boundaries are discussed 
as well. The models considered contain all the qualitative feature of quantized linearized cosmo- 
logical perturbations for a Randall- Sundrum universe having an arbitrary expansion history, with 
the sole exception of gravitational gauge invariance, which will be treated in a later paper. 
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I. INTRODUCTION 



In this paper we consider the problem of computing the cosmological perturbations 
for the one-brane Randall Sundrum model in which we inhabit a (3+l)-dimensional 
brane embedded in a pure (4+l)-dimensional anti de Sitter bulk spacetime. A Z2 symmetry 
about the brane is supposed and the perturbations are treated to linear order. The only 
bulk degrees of freedom are the (4+l)-dimensional bulk gravitons or gravity waves, which 
propagate in the bulk [l[ [| . These are emitted, absorbed, and reflected by the brane as 
well as reflected, or perhaps, more properly, diffracted by the bulk, making the problem 
of predicting the cosmological perturbations in a braneworld scenario much more complex 
than its (3 + l)-dimensional counterpart. 

A vast literature exists on the problem of calculating cosmological perturbations in a 
braneworld scenario. Some progress has been made, either by resorting to various approx- 
imations or by considering special spacetimes on the brane, such as eternal pure de Sitter 
space, for which additional symmetry can be exploited. A non-exhaustive sampling of the 
literature can be found in ref. Ql |i| and the references cited therein. In this paper we 
develop some techniques that could be employed to solve the problem completely, for a 
general expansion history and without resort to any approximations. The work reported 
lere, which is most closely related to considerations in Gorbunov, Rubakov and Sibiryakov 
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and the approach presented in the series of papers by Mukoyama [6], constitutes a 
first step toward this goal. Although we only discuss the connection with the Randall- 
Sundrum cosmologies, the methods and ideas developed here would apply equally well to 
more complex and perhaps more realistic situations with additional degrees of freedom in 
the bulk. 

It is this coupling between the brane and bulk degrees of freedom that renders the 
problem of cosmological perturbations in the braneworld scenario difficult. For linearized 
cosmological perturbations in the ordinary (3+l)-dimensional cosmology, the problem may 
be diagonalized by expanding the three spatial dimensions in Fourier components exp[ik-x] 
(where for simplicity here we assume a spatially flat universe) jjj. To linear order, there is 
no mixing between different k. Consequently, each such block may be analysed separately. 
The equations may be diagonalized further by separating the scalar, vector, and tensor 
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sectors. Here we define any quantity expressible in terms of derivatives acting on a scalar 
as scalar. Likewise, any quantity expressible as derivatives acting on a pure vector potential 
is regarded vector, etc. Under these definitions, spatial differentiation does not mix these 
three sectors, giving the evolution equations a block diagonal form. For each k there is a 
finite number of degrees of freedom whose time evolution is described by a finite number 
of ordinary coupled differential equations in cosmic time t. 

In the presence of an extra dimensional "bulk," the number of degrees of freedom in 
each k sector is greatly — in fact, infinitely — enlarged. (See, for example, the review 
When the size of the extra dimension is finite, the additional bulk degrees of freedom are 
discretely infinite, their spacing inversely proportional to the size of the extra dimension. 
When the extra dimension is very small, a large mass gap must be overcome in order 
to access the infinite tower of bulk excitations and at low energies the (3+l)-dimensional 
limit is obtained. However, when the extra dimension is infinite, the bulk degrees of 
freedom are continuous, labeled by an index k%, where < k§ < +00. For the case of a 
warped bulk spatial geometry, it has been shown how the coupling to the brane of the 
lowest energy degrees of freedom is suppressed, giving on large scales ordinary (3+1)- 
dimensional gravity on the brane. However, we would like to explore the corrections to 
such a limit. We would like to discover any possible differences in the predictions for 
the cosmological perturbations between the standard (3+l)-dimensional cosmology and a 
braneworld cosmology. Consequently, we must solve the interacting brane-bulk system. 

The vast enlargement of the number of degrees of freedom renders the problem difficult in 
two respects: Firstly, the equations become more complicated. But this is not all. Secondly, 
a solution to the problem of braneworld cosmological perturbations requires specifying 
initial conditions for the infinite number of extra degrees of freedom. While the first 
problem is of a technical, computational character, the second is of a more fundamental or 
physical nature. 

The problem of initial conditions for the bulk degrees of freedom can be further elu- 
cidated by considering the Penrose diagram for Randall-Sundrum Universes with various 
expansion histories, as indicated in Fig. Panel (a) indicates the patch on the conformal 
diagram covered by the standard Randall-Sundrum coordinates, having the line element 

ds 2 = dx\ + exp[-2x 5 /£] • \-dt 2 + dx\ + dx\ + dxf\ . (1) 
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FIG. 1: Various braneworld cosmologies embedded into the conformal diagram for 
maximally-extended AdS 5 . Panel (a) shows the triangular portion of maximally-extended 
AdS covered by the Randall-Sundrum coordinates. Time runs upward and the horizontal direction 
represents the fifth dimension, the other three transverse dimensions being orthogonal to the page. 
The family of lines focusing to the vertex on the right represents surfaces of constant Randall- 
Sundrum time. These asymptote in the future and in the past to the null surfaces H/ + \ and 

H(—\, which are the future and past Cauchy horizons of this triangular patch. The vertical curves 

connecting the lower vertex to the upper vertex of the triangular regions represent surfaces where 
the Randall-Sundrum coordinate for the fifth dimension is constant. These are also surfaces on 
which the scale factor a{x^) for the three transverse spatial dimensions is constant. As one passes 
from right to left this scale factor increases. In panel (b) the trajectory of the brane (indicated 
by the heavy timelike curve) for a static Randall-Sundrum universe, having a Minkowski induced 
geometry, is indicated. In panels (c) and (d) we consider expanding universes, where the motion 
of the brane to the left (with respect to the surfaces of constant X5) implies expansion. In 
panel (c) the brane worldline for a de Sitter induced geometry on the brane is illustrated. The 
brane emanates from the lower vertex, as for the static brane, but because of the rapidity of 
the expansion, strikes the boundary at conformal infinity before reaching the upper vertex. In 
other words, de Sitter proper time on the brane becomes infinite at finite Randall-Sundrum time. 
Panel (d) illustrates a braneworld universe that is initially inflating but whose expansion later 
slows down to become a dust-dominated universe. In this case, because of the deceleration, future 
infinity on the brane corresponds to infinite Randall-Sundrum time. In all cases there is a past 
Cauchy horizon Hr_\ on which initial data for the bulk modes must be specified to completely 

determine the subsequent evolution of the coupled brane-bulk system. 

Panel (b) indicates a static M 5 brane universe. Here the brane trajectory, indicated by 
the heavy curve, coincides with one of the x 5 constant surfaces. if(+) and are the 

bulk horizons with respect to the observers on the brane. Panel (c) indicates a Randall- 
Sundrum universe with a dS 5 geometry on the brane, and panel (d) indicates a universe that 
is initially inflationary but then reheats to become a decelerating (e.g., dust-dominated) 
universe. In all cases, there is the past horizon #(-), of particular interest to us, where 
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FIG. 2: Fundamental processes for the bulk-brane interaction. The thin vertical solid line 
indicates the position of the brane and the bulk is situated to its right. A bulk graviton, indicated 
by the dashed line, may strike the brane, annihilating itself and creating an excitation of the brane, 
indicated by the heavy vertical line, as illustrated in panel (a). Likewise, an excitation of the 
brane may decay, emitting a bulk graviton, as illustrated in panel (b). Because we are interested 
here only in the linearized theory, there are no vertices joining more than two external lines. 
Propagation in the bulk may involve scattering and is not in general along straight trajectories. 
Moreover, bulk gravitons may reflect off the boundary without changing the excitations on the 
brane. 

initial data for the bulk gravitons must be specified. 

An initial state for the brane-bulk system may be characterized completely by specifying 
the quantum state of the incoming gravitons on Ht\ and of the degrees of freedom on the 
brane at the intersection of the brane with Hi-\. Subsequently, the bulk and the brane 
degrees of freedom interact. The relevant fundamental processes are illustrated in Fig. |21 
As shown in panel (a), bulk gravitons may be absorbed and transformed into quanta on 
the brane. Similarly, as shown in panel (b), brane excitations may decay through the 
emission of bulk gravitons. These may either escape to future infinity, leading to a sort 
of dissipation, or be re-absorbed by the brane, leading to what from the four- dimensional 
point of view appears as nonlocality. Such dissipation and nonlocality are effects that cannot 
be incorporated into an effective four-dimensional effective field theory description, and it 
is hoped that by studying these phenomena one might be able to find some distinctive 
observational signature for the presence of extra dimensions. 

In this article we present some simple toy models where a single or finite number of 
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degrees of freedom is coupled to a one-dimensional continuum of degrees of freedom. Such 
(l+l)-dimensional field theories coupled to a boundary having its own dynamics are the 
analogue of a fixed k-sector of the brane-bulk system. The degrees of freedom on the 
boundary represent the degrees of freedom on the brane and the degrees of freedom of the 
one-dimensional continuum represent the degrees of freedom of the bulk. We consider a 
number of examples of increasing complexity that exhibit all the qualitative features of the 
brane-bulk system except those concerning gravitational gauge invariance and gauge fixing. 
These questions shall be treated in a future publication. 

In Section II we study the simplest such system. A harmonic oscillator is coupled to 
a stretched string so that excitations of the oscillator may decay, or dissipate, emitting 
waves on the string that propagate to infinity. If there are no incoming waves on the string, 
which is classically possible but not quantum mechanically, the motion of the oscillator is 
described exactly by the equation of motion for the dampened harmonic oscillator. Simi- 
larly, incoming waves may excite the oscillator. In the analogy to braneworld cosmology, 
the oscillator represents a degree of freedom on the brane and the string represents the 
bulk gravitons. We work out the quantum mechanical description of this system and in 
particular of its effective description in detail. We show how the annihilation and cre- 
ation operators localized at the end of the string, a osc (t) and aj sc (t), do not have definite 
frequency but rather are the superposition of modes of different frequency having a cer- 
tain characteristic width of order the classical decay rate 7. It is this spread in frequency 
that causes the commutator [a osc (t), a\ sc {t')] to decay in modulus when \t — t'\ becomes 
large compared to We also consider several oscillators coupled to each other and to 
a string and the case where the mass density of the string is not uniform so that there 
are reflections, leading to nonlocality in the effective description for the oscillator degree 
of freedom. All the systems in this section have a time-translation invariance, so that a 
Fourier decomposition in time can be employed, greatly simplifying the problem. 

In Section III we consider cases where there is no time translation invariance, in which 
the coupling evolves with time. While the examples of the previous section are analogous to 
the static Randall-Sundrum brane, where the brane geometry is that of (3+l)-dimensional 
Minkowski space, the examples of section III are more akin to the expanding braneworld 
universe, where the couplings of the various modes to the brane evolve with time. In the 
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time translation invariant case, at any given finite time, the quantum state of a mode 
depends only on the quantum state in the bulk at past infinity. However, with the time- 
dependent interactions, cases may be contemplated where the state on the brane depends 
on a linear combination of the initial state on the brane and that in the bulk. Some 
explicit examples are worked out where the linear canonical transformation — a sort of S 
matrix between the "in" and "out" states — is calculated. The matrix linking the "in" and 
"out" states for the bulk is akin to the calculation in Gorbunov, Rubakov and Sibiryakov 
[^] where the Bogoliubov transformation between the "in" and "out" states is calculated 
for a brane having a dS 4 geometry due to a pure cosmological constant and nothing else 
on the brane. 

In section IV we consider the complications that arise when the wave equation on the 
string is such that the general solution can no longer be decomposed into left and right 
movers. Any sort of mass term or non-uniformity in the speed of propagation will render 
such a decomposition impossible by causing disturbances emitted from the brane to scatter 
and propagate back onto the boundary. There can be many such multiple reflections. In 
the case of interest to us, the fact that the bulk geometry is AdS 5 rather than M 5 causes 
such scatterings, which might also be described as diffraction, particularly for wavelengths 
of order the AdS curvature scale. We develop a perturbative approach to treating such 
scattering. The final result indicates how the bulk interaction in the toy model can be 
expressed in terms of integral kernels. The approach is similar to that sketched in the 
papers of Mukoyama 

Finally, in section V we consider the case where the boundary follows an arbitrary given 
timelike trajectory x(t). (Here X5 is simply shortened to the x of our toy models.) We 
show how the relevant Green's functions may be constructed using virtual sources whose 
strengths are determined by solving a Volterra integral equation of the second type. These 
techniques are relevant to the brane-bulk problem because the propagation of the bulk 
gravitons is most trivial in a coordinate system where the bulk is static but the brane 
moves as a result of the expansion of the universe. In section VI we conclude with some 
closing comments. 
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FIG. 3: Simple model for the brane-bulk interaction. A simple harmonic oscillator con- 
sisting of a mass m bound to a spring of stiffness k = muj 2 is coupled to a stretched string of 
linear density \x and tension r. q(t) represents the vertical displacement of the oscillator, and 
the field u(t,x) represents the vertical displacement of the stretched string. Here the oscillator 
represents a degree of freedom localized on the brane and the modes on the string represent the 
bulk gravitons. As the mass oscillates, energy escapes into the bulk, leading to dissipation. If the 
linear density of the string is not uniform, some of these modes scatter back, leading to effects 
nonlocal in time. 



II. EXAMPLES WITH TIME TRANSLATION INVARIANCE 



A. Quantum Mechanics With Dissipation : A Single Mass Coupled to a String 

After the linearized perturbations of the brane-bulk system have been decomposed into 
their Fourier components with respect to three transverse spatial dimensions, the evolution 
equations decompose into independent blocks, each labeled by the three-dimensional wave 
number k. Each such block has the dynamics of the (1 + l)-dimensional field theory coupled 
to one or more degrees of freedom localized on the brane boundary. In this section, we 
study a simple mechanical model of such a block. We consider a simple harmonic oscillator 
attached to the end of a string (see Figure OJ). In the absence of coupling to the string, the 
equation of motion for the oscillator is 

mq(t) + mu 2 q{t) = F(t) (2) 

where the external driving force F(t) vanishes. The quantum states are simply the usual 
simple harmonic oscillator ladder of states of energy E n = Tiuj{n + ~). 

We now introduce dissipation in a 'conservative' way — that is, we couple the above 
oscillator to a system with an infinite number of degrees of freedom in such a way that a 
classically excited state of the oscillator progressively radiates its energy into the continuum 
and this energy never returns. (Some other approaches to treating dissipation within the 
framework of quantum mechanics may be found in refs. jl^.) For concreteness, consider a 
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string of linear density \i and tension r. For simplicity, we set // = r for the string, so that 

d 2 u d 2 u 

W~dx^ = ' ^ 

If we attach the string to the oscillator at x = 0, the driving force exerted by the string on 
the mass becomes 

F(«)=^| . (4) 

Classically, we solve this system by postulating as an initial condition the absence of 
any incoming (or left-moving) waves on the string. Using the ansatz 

u(t, x) = C exp[iu(x — £)] (5) 

and setting 7 = r/m, we obtain from (J2J) where q(t) = u(t,x = 0) the quadratic equation 

uj 2 — u) 2 = ijui, (6) 

having the solution 

«-±Z±^£. (7) 

This solution, of course, diverges at large positive x and large negative t; however, these 
divergences are unphysical, because in constructing a Green's function, we match to a 
vanishing solution below a line of (t — x) < t' where t' is the time corresponding to the 
singular source located at x = at the oscillator. Classically, from an effective point of 
view, we could dispense with the string and its degrees of freedom infinite in number and 
instead simply use the equation 
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q(t)+l<l(t)+u 2 q(t) =F ext {t) (8) 



where the force exerted by the string on the oscillator at x = 0, F string (t), has been placed 
on the left-hand side and is now the dissipative term rwyq^), and F ext (t) is included to 
denote an external force on the oscillator other than that exerted by the string. The arrow 
of time here was set by our initial condition that there be no incoming waves on the string 
in the past. If we had instead postulated no incoming waves in the future, the sign of the 
7<j term would be inverted. It is important to note that, from a classical point of view, 
eqn. (JHJ) is exact. There are no corrections. 
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Quantum mechanically, the initial condition above does not make sense. Vacuum fluc- 
tuations of the string require the presence of at least some left-moving waves. Otherwise, 
it would not be possible to satisfy the canonical commutation relations. Consequently, 
we must proceed by some other means. We start by postulating a vacuum state for the 
string-oscillator system. Note that when coupled to the string in the above way, the oscil- 
lator degree of freedom disappears. It has become part of the string, because we have set 
q(t) = u(t,x= 0). 1 

For a fixed time dependence proportional to e~ tuJt , we find the normal modes of the 
string with the ansatz 



u u (t,x) 

where u is given in (J7|). We have 



sin 



COX — 0(U) 



-iuit 



(u) 2 - u 2 Ju u (t, x = 0) = 7-rp(*, x = 0), 



or equivalently 



(oj 2 — UJ 2 ) HI) 



10J, 



dx 



-JU cos 



where 7 = (r/m). It follows that 

0H 



sin 



-70; 



"(7/2) 



(9) 



(10) 



'ID 



(12) 



/yV' + O^-w 2 ) 2 v /( 7 /2) 2 + ujf 

where the latter approximation holds when 7 <C u>, in other words when the quality factor 
of the oscillator is very large. 

The field on the string may be expanded as 

-°° du 



u(t, x) 



sin 



ux - (f)(u) a(u>) e + a{ujy e 



where 





[d(uj),a(uj'y] = 5{uj-uj') 



(13) 



(14) 



1 To understand how this occurs, it is useful to consider the string as a continuum limit of a series of 
small masses (Am) — fi(Al) connected by springs of separation (Al) attached to a harmonically bound 
mass of fixed magnitude M at the end point. In the continuum limit, the frequency of the mode where 
the large mass oscillates relative to its neighbors approaches infinity. It is always of energy of order the 
cutoff. 
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The normalization factor follows from the canonical commutation relations for the field 
u(t,x). It follows that the quantum operator corresponding to the oscillator degree of 
freedom localized at the end of the string at time t is given by 



du 



sin 



d(u) e~ iU + d(uY e +luJt 



"fU 



Define 



where 



^TTTWyU ^2^2 + ^2 _ ^2) 

p 1 ' 2 ^) = N 



d(u) e~ iujt + d(uV e +iujt 



■yu 



We may define 



Because 



N~ 



du 



^ 2 u 2 + (u 2 — u 2 ) 2 

JU 



it follows that at equal times 



7 2 cj 2 + (oo 2 — u 2 ) 2 
du p l l 2 {u) exp[— iut] d(u) 

+oo 

du p{u) = 1, 
a 0S c(t),al sc (t) 



(15) 



(16) 



(17) 



as for an ordinary harmonic oscillator. However, in the general case 



Qosc(^)) ^osci^ ) 



du p(u) exp[— iu(t — t')], 



(19) 
(20) 

(21) 



which for unequal times is of a modulus smaller than one. For \t — 1'\ <C 7 _1 this is to a 
good approximation a simple phase, namely exp[— iu(t — t')], just as one would obtain by 
ignoring the dampening and pretending that the oscillator is uncoupled with respect to the 
string. However, when \t — 1'\ becomes comparable to or greater than 7 _1 , one observes a 
sizable diminution in modulus. Unless the density p[u) has some delta functions spikes, the 
Riemann-Lebesgue lemma implies that, for times sufficiently well separated, the operators 
almost commute. 

Let us consider the special case of a uniform string of small tension, so that the approx- 
imation on the far right-hand side of eqn. ()12j) holds where 7 <C u. In this case 

1 (7/2) 



p(u) 



7l( 1 /2) 2 +(LU-LUy 



(22) 
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and with very little error we may replace the interval of integration [0, +00) with the 
doubly-infinite interval (—00, +00), so that 



a osc {t),a\ sc {t) 



+00 



duj 



(7/2) 



TTJ-oo ( 7 /2) 2 + (W - Uj) 2 

exp[-iiu(t - t')} ■ exp[-(7/2)|t - t'\\. 



exp[— iu(t — t')] 



(23) 



For broad resonances, there are corrections to this formula, but the general qualitative 
behavior persists. 



B. Multiple reflections and non-locality 

We may also contemplate more complicated systems for which the classical behavior of 
the oscillator is no longer local in time, as for example as indicated in eqn. (JH|) . but rather 
an equation of the form 



q{t)+'yq(t)+Q i q(t) = - 

m 



F ext (t)+ f dt'G(t-t') F ext (t') 



(24) 



where the kernel G(t — t') represents the amplitude for a wave on the string generated at 
time t' reflecting back to strike the oscillator subsequently at time t. Such nonlocal behavior 
will arise in the brane-bulk system because bulk gravitons initially emitted by the brane 
will be scattered back to the brane by the curved bulk geometry. 

For concreteness, consider the case where rather than being uniform, the linear density 
of the string abruptly increases at x = L, so that at that point the propagation speed 
suddenly drops from q to c r as one passes from the left to the right. Physically, a portion 
of the outgoing wave is transmitted and another is reflected back toward the oscillator. 
Across the junction from left to right, we have the following matching rules (which follows 
from continuity of the amplitude and its first derivative) 

(x-L) 



cos 



sin 



uj- 



UJ- 



Cl 

[x-L) 



Cl 





(x 


-L)~ 




COS 


UJ 










/Cr 




(x 


-£)] 




j sin 


u — 








c r 



(25) 



so that the solution for x < L satisfying the boundary condition at the oscillator 



sm[ujx/ci + (f)(uj)} — sia[u(x — L)/ci] cos[<P(uj)+ujL/ci]+cos[uj(x — L)/ci\ sin[<fi(uj) -\-uL/cj\, 

(26) 
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for x > L transforms into 



(c r /q) • sin[cj(x — L)/c r ] cos[(p (u) + uL/q] + cos[u(x — L)/c r ] sin[0(cj) + uL/q], (27) 

which has the overall amplitude 

A = \j(c r /ci) 2 cos 2 [0(cj) + uL/ci] + sin 2 [<f)(uj)+LuL/ci]. (28) 

Since the string is semi-infinite, it is this amplitude alone that determines the normalization 
of the modes on the string. The amplitude on the interval [L, +00) infinitely outweighs 
that on [0, L\. It follows that the spectral density takes the form 

— sin 2 [0(c<j)] 
P ^' = (c r /q) 2 cos 2 [0(cu) + uL/d] + sin 2 [0(cu) + uL/q] (29) 

where N is a normalization constant. Qualitatively, the density has the form of an ordinary 
resonance (from the numerator) masked by a function with spikes at </>(u) + uL/ci = nn of 
a certain width. The closer the reflection coefficient R — (q — c r )/(q + c r ) is to unity, the 
sharper the spikes. In the extreme limit c r — > 0, the reflection is total, so that the system 
is effectively no longer semi-infinite but rather a finite cavity, with a discrete spectrum. In 
this case the spikes have a 5-function character. The presence of these periodically spaced 
spikes may be understood as follows. The normalization of the modes of a given u is 
determined by the amplitude of the wave at x > L. However, it is only near a resonance 
that this wave has a significant penetration into the region x < L to the left of the junction. 
Consider the Fourier transform of 



f(T) = J2 R H exp[i5n]exp[-iLj(T -2nL)]exp 



-l\T-2nL/ Cl \ 



(30) 



which represents the effect on the commutator of (multiple) reflections, with the real reflec- 
tion coefficient R and phase shift 8 describing the effect of a single reflection. This Fourier 
transform is simply 

+00 

p(uj)=p Q (uj) Y R lnl exp[in5]exp[+2inuL/ci] (31) 

n=— 00 

where po(u) is the Fourier transform of exp[— iuT] exp — (7/2) |T| . We may evaluate 

+00 j p2 

Y R H exp[inS]exp[+2inujL/ci} = 7 ^32) 

„_oo " 1 + R 2 - R(exp[+i(2uL/ Cl + 5)] + exp[-i(2wL/c, + 5) A 
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Similarly, using the fact that 4>{oj) ~ 7r/2 in the neighborhood of the resonance, we may 
rewrite the denominator in eqn. (|2H|) as 



(c r /ci) 2 cos 2 [<p (u) + ujL/ci] + sm 2 [(ft (u>) + uL/ci] (c r /ci) 2 sin 2 [ci;L/Q] + cos 2 [ci;L/q] 



' (l + C 2/c 2 ) + (l-c2/c 2 )cos[2^L/ Q ] 

; f r i \ (33) 

1 + Acos[2cjL/q] v ; 

where A = (cf — c 2 )/(cf + c 2 ) = 2R/(1 + R 2 ). We observe that eqn. (}33|) has the form of 
eqn. (jHOJ) with phase shift 5 = ir. 



The above example illustrates how for models having a time translation symmetry non- 
local effects are encoded in the spectral density p(u>). Here the multiple reflections alter 
the spectral density Pq{uj) in the absence of reflections through multiplication by the mask 
function in eqn. (|33|) having the profile of a picket fence. After an integer number of re- 
flection times 2L the phases of the various pickets interfere constructively causing a peak 
in the commutator. 

We may also contemplate examples where the linear density of the string n{x) varies 
smoothly so that there is an amplitude for reflection everywhere and these reflections 
interfere with each other to give a rather complicated response kernel G(t — t'). The method 
of calculating this nonlocal kernel for the brane-bulk system shall be the subject of sections 
IV and V. 

C. Interpretation in terms of Hilbert space 

Let us consider the physical interpretation of this diminution in magnitude of the com- 
mutator. Our Hilbert space H str ing has the structure of a direct product of an infinite 
number of harmonic oscillator Hilbert spaces. There exists an infinite number of particle 
species. In an abstract sense, let /(^) be a normalized square-integrable complex function, 
that is 



Each such function specifies a type of particle whose annihilation operator is given by 



1 



1 



2 




(34) 




(35) 
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which together with its conjugate a J- generates a single-particle subspace Ti.f C lustring- We 
may also consider the Hilbert space of "particle species" (distinct from the above spaces) 
endowed with a natural inner product: 

(f\g) = / dcu /» g{u). (36) 

More precisely, particle species are labeled by complex rays in this Hilbert space: any 
two such complex-valued, square-integrable functions differing only by a complex phase 
represent the same kind of particle. 

Without coupling to the string, the type of particle coupled to the oscillator does not 
evolve with time. a(t) = exp[— iut]a(t = 0), and as noted a difference in phase does not 
correspond to a difference in particle type. However, in the case of coupling to the string, 
the oscillator annihilation operator, because it is a superposition of string modes of differing 
frequency, evolves in this ray space of different particle species. Any quantum operator 0(t) 
localized at the oscillator at the end of the string at a time t may be expressed as a sum of 
products of d osc (t) and a\ sc (t). Any operator that cannot be so constructed is not localized 
at the end of the string. 

Physically, the fact that for large t, a osc (t) almost commutes with a osc (t = 0) reflects the 
fact that the particles created by aj sc (t) escape from the neighborhood of the end of the 
string and propagate down the string, eventually leaving no trace detectable by an observer 
able only to operate on the string end. Once the particle has escaped from the end of the 
string, no operator localized at the string end is capable of altering its state. 

Quantum mechanically, from the point of view of a quantum observer only coupled to 
the mass at the end of the string, calculating probabilities that involve "dissipation" — that 
is, the escape of energy, particles, and information down the string — requires summing over 
all such final states. That is, rather than considering vacuum-to-vacuum amplitudes we 
must sum over the particles not seen that escape down the string. 

Let us first consider the simplest possible process, the creation of a single oscillator 
quantum at t — U and the outcome of an attempt to recover this particle at a subsequent 
time tf. One possible outcome is that the particle is detected at t — tf. The amplitude for 
this process is 

A = (0\a osc (t f ) aL(^)|0) = exp[-iu(t f - U)] exp[-( 7 /2)(t / - U)]. (37) 

15 



Note that the probability 

p lf ^ h = lAf^eM-litf-U)] (38) 

is strictly less than one, and 

Po f <-u = 1 - Pi f ^u = 1 - exp[-7(t/ - U)] (39) 

gives the probability that no particle is detected at the end of the string. This latter 
probability is actually the sum of the probabilities for an infinite number of orthogonal 
processes. One takes the infinite sum over / where the index / labels all the particle types 
orthogonal to the type coupled to the string at t = tf. In other words, 

p 0f ^ u =J2\A If <- U \ 2 . (40) 

More complicated 'dissipative' processes may be calculated analogously. 

III. EXAMPLES WITHOUT TIME TRANSLATION INVARIANCE 

In the previous section we considered harmonic oscillators coupled to a string. Our 
principal tool was the decomposition into Fourier modes of fixed frequency with respect 
to time. To quantize we associated the positive frequency modes with the annihilation 
operators of the unique preferred vacuum state. In this section we generalize to situations 
where a global time translation symmetry is lacking. We are in particular interested in 
cases where the coefficients characterizing the coupling of the oscillator or oscillators to the 
continuum evolve with time. In this case Fourier methods cannot be used. Instead it is 
necessary to work in real space. In a braneworld cosmology the expansion of the universe 
on the brane renders the system time dependent. 

For the most general case, it is not possible to single out a unique preferred vacuum. 
However for a massless wave equation, for which any solution may be decomposed into 
left-moving and right-moving waves, this is not a problem because it suffices to postulate 
the usual vacuum state for the incoming waves. In other words, where uj > the modes 
having the form exp[— iu(t + x)] in the limit t — > — oo are associated with annihilation 
operators and the exp[+icu(t + x)] modes with the creation operators. This defines a unique 
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vacuum state, even when the couplings to the harmonic oscillators at the end of the string 
and of the harmonic oscillators themselves vary with time. The same applies when the 
equation becomes massless sufficiently quickly as x — > oo. With respect to the oscillators, 
it is necessary to define an initial quantum state if the oscillators decouple from the string 
sufficiently quickly as t — > — oo. If they do not decouple rapidly enough, the initial state 
becomes erased because the initial state of the string degrees of freedom are infinitely more 
relevant. 

In the first worked example we consider a single oscillator whose spring constant, which 
is proportional to u 2 (t), varies with time. If ou 2 (t) approaches a constant as \t\ becomes 
large, there are well-defined "in" and "out" vacua. However, these vacua do not coincide 
with each other. The time variation in u) 2 (t) excites quanta through parametric resonance. 
The production of quanta may be characterized by a Bogoliubov (i.e., linear canonical) 
transformation relating the modes of positive and negative frequency for the "in" and the 
"out" vacua. In this worked example we consider the special case where uo 2 (t) has the 
form of a step function. In this case, the applicable integrals may be evaluated explicitly. 
However, the generalization to an arbitrary form for w 2 (t) approaching a constant for large 
\t\ is straightforward. 

The second worked example is analogous except that rather than one there are two 
harmonic oscillator degrees of freedom at the end of the string. In this example, for |i| > T 
the two harmonic oscillators are uncoupled. Consequently, to define an initial state, in 
addition to requiring the absence of incoming waves it is necessary to characterize the 
ground state of the oscillator initially uncoupled to the string. The resulting Bogoliubov 
transformation from "in" to "out" states acts on a Hilbert space with one discrete degree 
of freedom and a continuum of string degrees of freedom. The generalization to many 
oscillators or to more complicated time dependences is straightforward. 

A. A single harmonic oscillator with time varying spring constant 

We return to example of section IIA where the stiffness of the oscillator a) 2 , formerly 
constant, now becomes a time-dependent function u 2 {t). We consider the simplest system 
for which nontrivial Bogoliubov coefficients may be calculated explicitly. Concretely, con- 
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sider a harmonic oscillator of mass m whose coupling to a string of tension r is turned on 
for only a finite interval of time. The equation of motion is 

d 2 



m 



df^^ 



q(t) = F(t) 



(41) 



where the driving force F(t) exerted by the string on the oscillator is given by 

^/ \ du 



(42) 



x=0 



We may decompose the string solution, which obeys the massless wave equation, into the 
two components 



u(t, x) = u out (t -x) + u in (t + x). 



The solution for the oscillator is simply 



q(t) = u out (t) +u in (t). 



It follows that 



. . du 

~ T ~dx 



T 



'du in (t) du out {ty 



x=0 



dt 



dt 



(43) 



(44) 



(45) 



With 7 = (r/m), the equation of motion takes the form 

Uoutit) = (-) 

or more symbolically 



d d o / \ 



Uinit) 



Formally, 



D out {t) u out (t) = (-)D in (t) u in (t). 



U^t) = (") /* dt' G out {t | f') An(0 U in (t') 

J -co 



(46) 



(47) 



(48) 



where G out is the retarded form of D~^ t . This means that G out mediates the propagation 
from t' to t of the state u in subject to the potential D in . 

Consider the inverse Fourier transform of u out , which can be interpreted as the descrip- 
tion of the state u out in the Heisenberg representation. We obtain 

/ dtr , /'' 

111 J-oo 



— \ r+oo rtf f+oo 

u out (u f ) = — dt f e+^'f / dU G out (t f | U) D in (U) / dw, e~'^ u u in {u>i) 
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dui S(uf\ui) u in (ui) (49) 

J — oo 

where 

S(u f \u)i) = — / dt f e+^f / dii G out (tj | t,) An(*i) e-^^ (50) 

Z7T J — oo </ — oo 

is the corresponding propagator in momentum space describing the amplitude that a quan- 
tum with frequency u>i be later observed to have frequency ujf. If we associate e~ luJt with the 
annihilation operator, then S(uf\iVi) describes the probability amplitude that a quantum 
of frequency oji be destroyed and a quantum of frequency Uf be created. There are no 
homogeneous solutions for u out . u out is completely determined by u in . 
In order to work out a concrete example, we set 

j( 7/2)2 + ff , if | t |>T, 

\ ( 7 /2) 2 + n 2 , if |t| < r. 

It follows that the homogeneous equation D out (t) u(t) = has the solutions 

A exp[(- 7 /2 + iQ)t] + B exp[(- 7 /2 - itt)t], for \t\ > T, 



u(t) = < 



(52) 

C exp[(-7/2 + in)t] + D exp[(- 7 /2 - in)t], for |i| < T. 



We adopt a sort of Heisenberg representation where, here for \t\ > T, the coefficients of 
the two solutions are expressed with respect to an arbitrary origin at t — t' (or t = t")\ 

u (t) = A t , e ^/^m-t') + Bt , e (-7/2-<n)(t-f ) 

= A t „ e (-t/2+im-t") + Bt „ e (- 7 /2-in)(t-t"). (53) 

It follows that a change of origin is effected by the transformation 

(A\ (A v ,\ (A\ {e {-i/2+im"-t') o \ (A\ 

' ' !>(/" |0 

Similarly, for |£| < T, 



sy t „ V^w v ' \b) v \ o e (-7/2-ffl)(t"-f)j \ B ; t , 



(54) 



U (i) = C t , e (-7/2+i«)(t-f ) + A , e (-7/2-^)(t-t') (55) 

and a change of origin is effected by the transformation 

(C\ fC\ / e (-7/2+^)(t"-f) o \ /C\ 

Ul = P(t *^U),= ( o e.-^-nJlJ,/ W 
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Matching at a jump in u 2 , as (Q <— Q), is effected by the transformation: 

(C) =J M) = W<? + "> _ (57) 

V£>A +e UA_ e 2^V(fi-fi) (0 + ^)7 \Bj t _ e 

and similarly, in the other direction, as (Q <— Cl), 

( A ) =J-i( G ) = J_f( n + ") (5g) 

In order to construct Green's functions, it is useful to divide the the real line into three 
regions 

region I = (— oo, — T], 
region II = (— T, +T), 

region III = [+T,+oo). (59) 

A Green's function for the oscillator where both endpoints lie within region I (i.e., 
U < tf < — T), is expressed in terms of the conventional notation as 

G(t,U) = (d' + ^ + uj 2 )- 1 ^) = Sm[Q( ^~^ )] exp[-( 7 /2)(t-^)] 9(t-U). (60) 

In terms of our notation, where the homogeneous solutions are expressed according to the 
conventions of eqns. (J53|) and (J54j) . the Green's function is generated by the column vector 



+i 

-1 



(61) 



localized at £ = U and placed at the far right in matrix expressions time-ordered from right 
to left. It follows that, in terms of our matrix notation, the Green's function is re-expressed 

as 

G(t f \t l )=( 1 ) P(t f \t t ) (62) 
^ 1 ' ^ 2in ' 

For propagation from region I into region II one instead would have 

G(t f \U) = (]) P(t f | - T) J P(-T \ tl ) (63) 
^ 1 ' ^ 2in ' 

and the other cases may be worked out analogously in a straightforward manner. If ti lies 
in regions I or III, the column vector (J61)) is used; otherwise, if ti lies in region II, Q in eqn. 
(jfiTjl is replaced with Q. 
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We now proceed to express the S matrix as the following sum 



S(uj f \uji) = Si<-i(uf\iVi) + S n <-i(uf\u;i) + S n i<-i(vf\ui) 

+ S n <-n(vf\vi) + Sni<-ii(vf\vi) + Sin<-iii(vf\vi), 



whose six terms are given by: 



S M (w f \u;i) = — (-) 



[ T dt f [ tf dU e +^ftf e -^ti 

<J — OO J —CO 



X 



+1 
2iU 



P(tf\U)[2 
■ ' ^ 2in 



1 

27' 



r 



dti 



-T 



-T 



T 



x J P(t f \-T) JP(-T\U) 



+i 

2ifi 
2ifi 



Siii^i(u f \ui) = 



-j- + fi 2 - w 2 + «7Wj 



J+T J-co 



2«Q 

-1 
2ifi 



>II*-II 



(Uf\Ui) = — (-) 



7' 



+ fr - a; 2 + i7^j 



J T dtf J dU e+^&e 



■iujit. 



x Pi/, /;)' 



-1 



S , 777 < -7/(w/|w i ) = 



-j- + fi 2 - ^ 2 + «7^j 



+T 



x ( J P(^ |+T) J" 1 P(+T | i 



+oo /-+T 

dtf J dU e+^&e-™' 

( 2itt \ 



S ni +-ni(wf\wi) = ^(-) 



r 



+ il 2 - uo 2 + i^uji 



/ dtf / d^ e +«*/t/ e -«* 

J+T J+T 



1 



xii P(i / 1*0 



+1 

2iQ 

2if7 
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We next proceed to calculate these six terms more explicitly. First, for Sj<-i we may extract 
and evaluate the integral 

/ dt f / dU e +iuj ^ e~^ iU 

J — oo J — oo V 



e (-7/2-in)(t/-*i) 



f-T , , ,+oo /e (-7/2+in)r o 

= / dt f e +*(«/-«i)*/ / dr 

J-oo JO 



e (-7/2-«i)T y 

W %> V U (7/2+^-1^) ' 

For Sjj^j, we may evaluate the two integrals: 

df e -^U e (-7/2±^)(-T- ti ) = e +^T 1 (72) 

J-oo 7/2 =F ill - iuji 

and 

/ dt f e^f e (-^m f +T) = £ e _ — J_ (73) 

i-T 7 7/2 ^ifi-iw/ 

For Sju^ii we have the first integral from again and also 

/ <ft f e (-7/2±«>)(t/-T) = e +-/T f (74) 

Sii<-u has the two nested integrals 

/■+T /-t, 



1 dt, y ' eft* e (-7/2±in)(t/-*i) 

= ^ (it/ e (-7/2±in+io; / )t / ^. e (+ 7 /2=Fin-jo; i )t j 

/•+T - p (j/2^in-iu>i)tf _ p -(-Y/2=fi£l-iu>i)T 

= / tft/ e (-7/2±«l-H«/)t/ 

1 



7/2 =1= ifl — icu, 
Siii<-ii has the integral 



e +i(u f -Wi)T _ e -i(w f -Wi)T e [2(-7/2±in)+i(a; i +a; / )]T _ e +i(u>i-u f )T 



(.75) 



hT - p-i^T _ +iu>iT (- 7 /2±ifi)2T 

dii e"^ e (-7/2±^)(^+T) = £ £ ^ (76) 

-t 7/2 =f ifi - iu;. 



Sjjj^in has the integral 



dt f / dii e + ^*' (r^ u e (-7/2±<n)(t/-*i) 

-T V+T 

/+00 rtf—T 

dtf e+^f-^f J dr e +*"iT e (-7/2±in)T 
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We have used the fact that 

r+oo 



Jo 



dt e +iuJt = 



+i 



uj + it 

For further reference, we note that 

1 1 



o 



— i 



dt e +iuJt = 

oo UJ — It 



(78) 



2tt8(uj) = 



y UJ + it uj — It 

where P denotes principal part. 
The final S matrix takes the form 



P U) 2 



1 1 

+ 



uj + it uj — it 



(79) 



S(ujf\uJi 



-28 (ujf — uJi 



nonsmg 



(uj f \uji 



}-n 2 - 7 2 /4 - i-yuf 
ujj — Q 2 — 7 2 /4 + i^ujf 



(80) 



where the singular 5-function term reflects what S would be if Q = Q (i.e., if the string 
constant of the oscillator were unchanged) and S nons in g is the smooth part resulting from 
the Fourier transform of the wave packet emitted as a result of the change during the 
interval [-T,+T\. 

The non-singular contribution to the final S matrix is given by 

~ 'uj 2 - Q 2 - 7 2 /4 - i^uJi\ ( ujf - Q 2 - 7 2 /4 - ryc^\] sin [(a;/ - uJi)T] 



s uj 2 -n 2 -7 2 /4 + i7Wi 



+ 



uj 2 — Vt 2 — 7 2 /4 + i^uji 

7^(Q 2 - n 2 ) 



7l(ujf - UJi) 



IT 



X 



{{Ui + ? 7 /2) 2 - W}{(ui + ii /2) 2 - n 2 }{(uj f + * 7 /2) 2 - n 2 }{(uJ f + ii/2) 2 - Q 2 } 



(uj f + uJi + i'yj x j- 



exp[+i(uj f - uJi)(-T)\ 



uJj + tl/2) 2 

+ [(uj f + i-f/2) 2 - tt 2 \ exp[+i(u f - uJi)(+T)} 
+ exp[- 7 T] exp[+iuj f (+T)] exp[-^(-T)] (Q 2 - Q 2 ) 

x^ujf + uji + i^cosp^T)} -i (uj i + i'j/2)(ujf + ii/2) + tt 2 



sin[ft(2T)] | 



n 



.(81) 



The term on the first line arises from the fact that the reflection from the oscillator during 
the interval [— T, +T] occurs with a different phase from that during its complement when 
\t\ > T. This term, which approaches a delta function as T — > oo, assumes that during the 
interval [— T, +T] the oscillator is described by its asymptotic amplitude and phase (as if the 



23 



change from Q to Q had taken place in the infinite past and there were no transients). This 
picture, of course, is an approximation because there were transients, and the remaining 
terms give the form of these transients. The two terms on the third and fourth lines have 
the form of scatterers localized in time at t = — T and t = +T, respectively. The first is 
due to the transient as the oscillator changes its amplitude and phase just after t = —T 
as a result of the change from f2 to Cl. The second term is the same at t = +T as the 
oscillator frequency changes back from Q to Q under the assumption that the previous 
initial transient has completely decayed away. Finally, the terms on the last two lines 
result from the fact that the first transient from t = —T has not completely decayed away 
by the instant t = +T. 

The transformations just calculated may be cast into a more general framework. What 
we have just been doing was calculating the matrix elements of a linear symplectic trans- 
formation relating the description of the system in terms of "in" modes to an equivalent 
description in terms of "out" modes. The two descriptions are related by a linear transfor- 
mation, which we may express more abstractly as 

fa out {uj)\ roo f S+ + {u)\u') S + -{u\u)')\ / a in (u')\ 

{al t (uj))~'° " U-+M"') S^(u\uj')) Ul n (u/)J' 
We know that both the "in" and the "out" operators must satisfy the canonical commuta- 
tion relations. In particular, for the "in" operators, we have 

[a in {u),a\ n {u')\ = 5(u-u/), 

[a in (uj),ai n (u')] = 0, (83) 
where u, u' > 0, and similarly for the "out" operators 

[a out (u),al ut (u')] = 8(uj-uj'), 

[a out (uj),a ou t(uj')] = 0. (84) 

However, the linear transformation (|82|) allows us to calculate the "out" commutation 
relation using eqn. (|83j). The requirement that the commutation relations thus obtained 
agree with those in ()84j) leads to a number of consistency conditions for the coefficients in 
(|%2*j). which we give below. 

By taking Hermitian conjugates, we obtain that S a p(uj\uj') satisfies the conditions: 

S ++ {uj\uj r ) = SI _(u\ui'), 
24 



S+_(u;|u/) = S*_ + (co\u'). (85) 

Eqn. (jHl|) gives 

/•OO _ , 

/ duo" [S ++ {uo\uo")Sl + {uo'\uo") - S + -{uj\u")SI_{u'\u")\ = 5(oj - w') (86) 

and 

/•oo 

dw" [S ++ {u\uj")S + 4uj'\uj") - S + 4u\u")S ++ {u'\u")) = 0. (87) 







These conditions are necessary and sufficient. 

For the "out" states or quanta, we may construct a number density operator N out (u>) = 
<^out{u)ciout{u) indicating the density of the number of particles of frequency u in the "out" 
space. We find that the expectation value of N out (uj) for the "in" vacuum state is 

(0 in \N out (u;)\0 m ) = ^'|5 + _(^|^)| 2 - (88) 
J o 

It follows that the "in" and "out" vacua will be different as long as 

/•OO /"OO /"OO 

(N out )= duj (N out (uj)) = duj dLu'\S+4uj\uj')\ 2 > 0, (89) 

JO JO JO 

that is, as long as the expectation value of the operator N out = J duoN out {uj) for the total 
number of "out" particles created from the "in" vacuum state is different from zero. We 
will obtain a vanishing expectation value, hence coinciding vacua, in the case when the 
string-oscillator coupling remains constant. 



B. A pair of harmonic oscillators with time dependent parameters 

In this section, generalizing the techniques of the example, we solve the equations for a 
pair of coupled harmonic oscillators with a forcing term included 



m 



dt 2 + 



'cu 2 (t) a 2 it) 
.a 2 it) oj 2 (t) 



(90) 



Here a 2 (t) is the coupling between the oscillators, which we shall make approach to zero 
sufficiently fast as \t\ — > 0. We attach only oscillator 2 to the string so that 



qi(t) = q OS c(t), q 2 (t) = u out (t) + u in (t) 



(91) 



25 



where the solution on the string is 



ll(t, x) = U out {t -x)+ U in (t + X) 



for x > 0. It follows that 

*U*) = 0, F 2 (t) = rf x 
The equation of motion takes the form 



x=0 



[dt 2 + {q 1 ) dt + { a 2( t ) u 2 



D\t) 
D 2 {t). 



U in (t) 



where we set 



-£ + -ri-«t. 



a 2 , if \t\ < T, 



'du in (t) _ du out (t) \ 
dt dt ) 

( q sc{t) 
\u out (t) 



-a 2 (t) 



0, if |i| > T. 

For \t\ > T, we may write the homogeneous solution in the form 



where 



4o e ~ 1(a)t 



and 



+ Z7(i) = W(i) = 

+17(2) = ^(2) = -<^i, 

+*7(3) = w (3 ) = +i[7/2 + v/(7/2) 2 -^ 2 2 



+«7(4) = ^(4) = +*b/2 - \/ (7/2) 



2 " ^ 2 ] 



^(1) - ^(2) 



1 







or ^) = ^) = li 



For |i| < T, the u>(s,)S are the roots a; of 



—a 



—a 2 w 2 — a;| + rycj 



0, 



or equivalently of the quartic polynomial 

(to 2 — oj\)(uj 2 — uj\ + 270;) — a 4 = 0. 
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[Note that for 7 7^, a 7^ this quartic does not reduce to a quadratic in u 2 .} The corre- 
sponding eigenvectors are given by Efy, so that 

E{a) e-" (a)t (101) 

form a complete basis for the homogeneous solutions. The indices (a) refer to the basis for 
\t\ < T whereas the indices (a) refer to the basis for \t\ > T. 
We define the matching matrix J(a)( a ) so that for all (a) and i 

E (a) = Y, J (a)(a) E (a), 
_ (a) 

%) E (a) = J2 J (a)(a) l(a)E\ a y (102) 
(a) 

We define G 1 ,^ such that for all i, j 

E G\ a) E{ a) = 0, 

(a) 

T,7(a)G\ a) Ei a) =5 ij: (103) 

(a) 

and the G\ a \ are analogously defined. 

Consequently, for t, t' entirely in region I (or entirely in region III) the Green's function 
is given by 

G^{tX) = E\ a) P {a){b) {tX)G\ b) , (104) 

where 

P (am (t,t>) = 5 {am e-TwC"*') = P {a) (t,t') S (am (105) 
describes a change of origin, and similarly for points entirely in region II, 

G ij (t,t')^El &) P {h){i) (t,t')G^ y (106) 
When points straddle the regions, example for propagation from / into II one obtains 

G*"(f,f ) = Efa P (a)(5) (f, —T) J m P ma) (-T,t') G\ a) . (107) 
We now give the expressions for the various Bogoliubov coefficients. 

^(OSC, +,OUt <- OSC, +,in) = ( J _1 )( a=2 )(a) P(a)(+T,-T) J(a)(b=2) (108) 
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S(osc, out <- osc, +, in) = ( J -1 )^^) P ( a)(+T, — T) J ( a)(&=2) (109) 

S(oSC, +, OUt <- OSC, -, in) = ( J _1 )( a=2 )(a) A&)(+ T > ~ T ) J («)(6=l) ( 110 ) 
^(OSC, -,OUt <-OSC,-,in) = (J _1 ) (a=1)(a) P ( a)(+T,-T) J(a)(6=l) (HI) 



S'(osc, +, out <— a;/, in) 

-T 

— ia;/t/ 



/ dtj (J- 1 )(a=2)(a) P(a)(T, —T) J (a)(6) P(ft)(— T, t/) GP {b) rP {tj) e 
J—oo 

+ eft, ( J- 1 )(a=2) (a ) P(«)(+T, f,) Gfo^fa) e -*"" (112) 



«S(osc, — , out <— a;/, in) = 

r T d*/ (^ 1 )(a=l)(a) P { a)(+T,-T) J (a){6) P( 6 )(— T, t/) G^P^j) e"*"'' 
+ l^rft, (^ 1 )(a=l)(a) P { a)(+ T ^l) {tj) (113) 



S(u>f, out <— a;/, in) = 

/■+oo /■- T 

/ tft F / df , Pg 2 P (a) (t F , +T) (J- 1 ) (a)((i ) P(a)(+T, -T) 

J -\-T J—oo 

S) 



xJ ( a)(6) P(b)(—T, tj) G} b) D*(tr) e-*"*' 



r+oo r +T 

+ <**f y r d* 7 P ( t } 2 P (a) (t F ,+T) (J- 1 )^)^) P(a)(+T,t/) Gja)^^/) e-*"«' 

/+oo r+oo 

dt F e +i "^ J t dt / P ( t ) 2 P (a) (t F ,+T)Gj a) P^(t / )e-^ 



/+T J+T 

+ e + ^*- y^dt/ Pgf P { a)(t F ,-T) J (a)(6) P ( 6)(-T,t 7 ) G j ^D j (ti) e 

-iuiiti 



-iojjtj 



+ T T dt F f ^ Pg 2 P(a)(tF,t/) GL^^j) e^'*' (114) 

J—oo J—oo 



S(u F , out <— osc, +, in) = 

/ dtp e+*"* Pg 2 P (a) (t F , +T) (J- 1 )^) P(o)(+T, -T) J (a)(6=2) 

+ /^f e +iw ^ Pg 2 P (a) (* F , — T) J (a)(6=2) (115) 



S(ujp, out <— osc, — , in) = 

/+oo 
^ di F e +i " FtF Pg 2 P(a)(tF,+T) (J" 1 )^) P(a)(+T,-P) J(a)(6=l) 
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+ / dt F e +t ^F El=2 p ( . )(tF; _ T) j (a)(fe=1) . (116) 

For the string mode, a positive frequency for uj,u>f corresponds to the annihilation op- 
erators, a negative frequency corresponds to creation operators. In the above formulae, 
summation over repeated indices is implied, and in some cases indices are repeated three 
times because intervening Kronecker deltas have been suppressed. 

The Bogoliubov coefficients connecting the oscillator to itself, as given by eqns. (jl(J8|) - 
(jlllj) . can be written as 

SfoSC, Sout *- OSC, S in ) = [Ss ou t,+$2,(a) + ^ out ,-^l,(a)] [Ss in ,+82,(b) + ,-$l,(6)] 

(J-Xxc) e"^ (2T) J( 5)(6 ) (117) 

where Si n , S out = +, — . For large T, the sum over (c) is dominated by the mode that decays 
most slowly among the eigenvalues of eqn. (jlOOjl . which are in general non-degenerate. We 
shall denote this mode as (cd) and its eigenvalue as T, so that in this case eqn. (|117j) above 
may be approximated as 

[5s out ,+52,(a) + S Sout ,-Si,(a)) [Ss in ,+^2,(b) + 6 Sin ,-6i,(b)) (^ _1 )(a)(c ) J(c D )(b) e~ T{2T \ (118) 

which approaches zero exponentially fast for large T. This attenuation can be interpreted as 
the decay of the oscillator excitations into the bulk. When (TT) ^> 1, any initial quantum 
information about the initial state of the oscillator becomes almost completely dissipated 
into the bulk. An initial state expressible as a tensor product of oscillator and bulk states 
becomes highly entangled. This means that if today we measure the state of the oscillator by 
recourse to operators acting only on the oscillator in terms of the initial state, we are almost 
exclusively measuring the correlations of the bulk modes — that linear combination that 
subsequently scatters and excites the oscillator. If the measurement is instead performed 
just before decoupling, the above holds with the straightforward modification that the 
second J matrix is suppressed. 

Oscillator excitations caused by incoming bulk excitations are described by the Bogoli- 
ubov coefficients in eqns. ()112j) - (jllHj) . We can rewrite these coefficients as 

^(osc, Sout <- in) 

= [iuA(.)+«uA(.)l (J^W) e- 7(£)(2T) J(S)(b) [^dtj e -Tw(" T -*') G^D^) e 
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+ [5s out , + 6 Ua) + 5 SouU -5 Ua) ] (J-% m " dtl e~ 7(E)(T_tl) G^D^tj) e 

J — T 

since D^ =1 (t) = —a(t) = for t < —T. Computing the integrals we obtain 



iujjtj 



(119) 



S{o^S out «- W/,in) = [<5so««,+*2,(o) + <*S out -<*l,(a)] (■/ ^(ojfg) 
/ e -7(E)(2T) giwjT _ _ g-iw/T _ e -7(c)(2T) e ia;rT _ 

I ; .7(g)(6) G^D 1 ' 2 + : 



G (a)^ 



(120) 



\ 7(6) - iUJ i 7(c) - ^/ 

For large T, we can suppress the exponentially decaying terms, which reflect the fact that 

the oscillator modes have not completely decayed into bulk excitations. We then find that 

eqn. ()120j) can be approximated for large T as 



[8s out ,+S2,(a) + $S out ,-Sl,(a)] {J %)(£) G^D 3 



-ILOjT 



(121) 



7(£) - i^i 

The Bogoliubov coefficients describing production of outgoing bulk excitations generated 
by oscillator excitations, eqns. (jll5j) - (jll6)) . also consist of two contributions: one from —T 
to +T during coupling and another transient after +T. We can rewrite as 



S(cop, out <— OSC, Si. 



dtp e 



iujptp 



+T 



+ / dt F e 46 "*' e-^ +T m^ J 5>6 [S Sin , + S Ub) + 5 Sin> _5 h 



(6)1 



e -7(c)(2T) e ia; F T _ _ e -ico F T _ g -7(c) (2T) e ioj F T v 



^(a) 2 ^(a)® — — +#l~ 2 

7(a) - 



'(c) 



7(c) - 



Xj(c)(6) [<fe, B ,+<*2,(&) +<fe 4 „,-A(6)]- ( 122 ) 

For very large T the terms suppressed by an exponential factor e~ 7 ( £ ) < - 2T - ) reflect a correction 
due to the fact that the oscillator mode has not completely been converted into bulk modes. 
With this correction suppressed, we find for very large T the following approximation for 
eqn. (fT22l 



-iw F T 

(123) 

7a - tto p 

which can be interpreted as a resonance of the oscillator for the bulk excitation produced. 

In the general case (compare with the end of the previous section), the linear transfor- 
mation between the "in" and the "out" operators may be expressed as 

/ d out {osc)\ /S' ++ (osc|osc) ^-(—(osclosc) S ++ (osc|u/) S+^oscIa/) \ / a lTl (osc) \ 



aL(osc) 



S_+(osc|osc) S (osc|osc) 5_ + (osc|u;') S (osc|a; r 

5 ++ (a;|osc) S , + _(cj|osc) S ++ (u\u') S+-(u\u') 
V 5_ + (a;|osc) S (c<j|osc) S-.+ (u\u/) S (w|tt/) J 



ai n {LU 
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(124) 



where summation and integration over repeated indices or variables is implied. In the case 
of several discrete oscillators, (osc) is replaced with an index taking integer values. 
The transformation must preserve the commutation relations 



[a iri (osc),aJ n (osc)] = 1, 

[a in (u)),a\ n (u))\ = 5(uj-uj'), 



and 



;o«t (osc), aj^ (osc)] = 
[a out (uj), ~ 



which imply that 



s ++ 


'osc osc) 


= s*_ 


_ (osc osc 


S+- 


^osc osc) 


= si 


+ (osc osc 




_(osc \uj') 


= S*_ 


_(osc \uj') 


S+- 


_(osc|u/) 


= S*_ 


+ (osc \uj') 


S ++ (u\osc) 


= S*_ 


_{uj\osc) 


s+ 


_(cj osc) 


= S*_ 


+ (uo osc) 


s. 


t+(w\uj') 


= S*_ 


_ (cj JcU ) 


s. 




= S*_ 





and 



(125) 



(126) 



(127) 



S ++ (osc | osc) S+ + (osc | osc) — S' + _(osc|osc)S'^_(osc|osc) 

+ J du" [S ++ (osc\uj")S* ++ (osc\uj") - S ++ (osc\lu")S* ++ (osc\uj")} = 1, (128) 

S' ++ (cj|osc)S'^ + (a; / |osc) — S , + _(o;|osc)S'^_(a; / |osc) 

+ J du" [S ++ (u\u")Sl4u'\u") - S+4u\u'')SI_(lu'\u")\ = 6(u-u'). (129) 

The remaining commutation relations 

[a iri (osc),a in (osc)] = 0, [a\ n (osc), d\ n (osc)] = 0, (130) 
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[a in (u), a in (uj)] = 0, [a] n (uj),al n (uj)] = 0, (131) 

imply that 

[a out (osc),a out (osc)] = 0, [a\ ut (osc), a Lt(osc)] = 0, (132) 

[a out (u;), aw(cj)] = 0, [aL(^), aLtM] = °> ( 133 ) 
which entail the additional conditions 

S' ++ (osc|osc)S' + _(osc|osc) — S' + _(osc|osc)S+ + (osc|osc) 

roc 

+ / du" [S ++ {osc\u'')S + -{osc\u'') - S + -{osc\u'')S ++ {osc\u'')} = 0, (134) 

JO ' 

5 ++ (u;|osc))S' + _(u/|osc) — S' + _(u;|osc),S' ++ (u/|osc) 

roo 

+ / duo" \S ++ {u\u")S + Ju'\u") - S+Ju\u")S ++ {J\u")\ = 5{u - J). (135) 

JO 

Similarly, we can construct the number operators for quanta in the "out" space, 
namely N out (osc) = ao Ut (osc)a OU i(osc) for the number of oscillator quanta and N out (u) = 
olut( U} )Q'<mt(w) for the number of bulk quanta of frequency to. We find the corresponding 
expectation values 

roo 

(0 in | (osc) 1 in ) = |S' + _(osc|osc)| 2 + / du'\S + _(osc\uj')\ 2 (136) 

Jo 

<0 in |i\U(w)|0 in ) = |S + _(u;|osc)| 2 + / du'\S+-(u\u')\ 2 , (137) 

Jo 

which will vanish if the "in" and "out" vacua coincide. 

In the application to braneworld cosmology, when we observe the cosmological per- 
turbations today, we will (using the assumption of Gaussianity) be measuring, in effect, 
expectation values of observables quadratic in the creation and annihilation operators lo- 
calized on the brane today, in other words, the "out" oscillator operators. Ignoring the 
additional discrete indices (which we suppress), we offer a useful parameterization of the 
relevant part of the S matrix. Let ab rane ^ out and al rane out be the operators of which we 
want to take the quadratic expectation value. The S-matrix expresses these as a linear 
combination of a 6rane;0Ui and a\ raneout , on the one hand, and of a 6ui/Mn (u;) and a\ ulkM (uj), 
on the other. The following offers a useful parameterization of this transformation. Require 
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that A(, rane j n and Abuik^n, normalized SUch that [^&rane,in> ^.Lane,m] ~ \A-bulk,im A\mlk,ini ~ ^' 

be entirely on the brane and in the bulk, respectively. Then ab ra ne,out may be expressed in 
terms of these according to one of the three following possibilities: either 

a brane,out = COS # A.brane,in + SU1 # ^bulk,in (138) 

where < 9 < tt/2; or , 

Q>brane,out = COsh £ Ab raney i n + sinll £ Al u i k ^ n (139) 

where < £ < +00; or 

^brane,out 

sinh £ Al aneM + cosh £ A bulkM (140) 

where < £ < +00. We observe that the bulk initial state may have a very important, even 
dominant, role in determining what we see on the brane today. Ab ranetin , and its conjugate, 
may in turn be related to a branetin by a Bogoliubov transformation far from the identity, 
and the same applies to Ab u ik,in- 

IV. GENERALIZATION TO NONTRIVIAL DYNAMICS ON THE STRING: RE- 
FLECTION, DIFFRACTION, AND DISPERSION 

In the previous section, because the undulations on the string were described by a 
massless wave equation, it was possible to separate u(t,x) on the string into a sum of the 
form u in (t + x) + u out (t — x). The left-moving component u in (t) acted as a source for the 
oscillator, which in turn radiated exclusively into the right-moving u out channel. However, 
for a massive wave equation of the form 

d 2 d 2 



W 2 ~dx 2 - +m2{x) 



u(t,x) = 0, (141) 



it is no longer possible to separate into left-moving and right-moving components. If we 
consider the m 2 (x)u(t, x) term as a perturbation, we see that initially right-moving out- 
going waves are scattered back onto the oscillator situated on the boundary, and incident 
left-moving waves, rather than striking the oscillator, are sometimes rescattered back to- 
ward infinity by this term. In the braneworld problem, the curved geometry, due to the 
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fact that the bulk spacetime is AdS or AdS-like rather than Minkowski, introduces such 
scattering. For the time-independent problem, all these effects can be treated straightfor- 
wardly by decomposing into Fourier modes of definite frequency. Nothing more is required 
than solving eqn. (j!41|) by separation of variables with the appropriate boundary condi- 
tions, following the formalism developed in section II. However, for the time-dependent 
problem, when no such decomposition is possible, it is necessary to solve for these multiple 
reflections in real space rather than in momentum space. One possibility would be to treat 
the m 2 (x)u(t,x) term as a perturbation that generates an infinite series expansion in the 
number of reflections. However, for all but some very special cases, this series is either 
divergent or very slowly convergent. 

A well-behaved series may be obtained by noting that for all frequencies except those 
near the resonance uo « uq, there is very little interaction between the string and the 
oscillator. Setting 



u ou t(u) = P(u) u in (u), (142) 



we observe that 



= (-) • u \-f-^ - d«) 

ur — uj z + z^lu 

Away from the resonance, the end of the string acts almost exactly as a Dirichlet boundary 
condition, for which one would have P{oS) = —1 exactly. Consequently, to obtain a rapidly 
convergent series, it is advantageous as a first approximation to solve for the propagation as 
if there were Dirichlet boundary conditions. Then the violation of the boundary condition 
for the oscillator becomes a source from which a perturbative expansion by successive 
approximations can be generated. 

More explicitly, for a given incoming wave, with Cauchy initial data specified on past 
null infinity J\f(~\ we first find a continuation Uo(t, x) such that the Dirichlet boundary 
condition 

u {t,x = 0) = (144) 
is satisfied. We shall call the Uo(t,x) obtained in this way the "incident wave." We then 



construct a series expansion 



u(t,x) = ^2u n (t,x) (145) 

n=0 
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such that 

(146) 

x=0 

is satisfied. Here the string-oscillator coupling 7 (equal to r/m where r is the string tension) 
serves as the expansion parameter. We correct for the violation of eqn. ()146j) by uo(t, x = 0) 
by choosing boundary data for U\ such that 

(147) 

2 = 

by means of a retarded Green's function for the oscillator. Then this Dirichlet boundary 
data for u\ on the boundary is used to extend U\ to the full domain x > using the 
Dirichlet form of the retarded Green's function in the bulk, in other words 

Ul (t, x) = / dt'—G D (t, x] t', x' = 0) ux{t', x' = 0). (148) 

J —oo ox 

This process repeats itself, so that the relations 

(149) 
and 

u n+ i(t, x) = f dt'^G D (t } x] t', x' = 0) u n+1 (t', x' = 0) (150) 

J —00 ox 

are satisfied. At the boundary we use the Green's function for the oscillator to obtain 

/t Q11 
dt' G osc (t, t') 1 -^{t',x' = (151) 
-00 ox 

In terms of propagators, we may write 
u(t,x) = u (t,x) 

+ 7 f dt' f dt® D x ,G bulk (t,x;t',x' = 0) G osc (t';t (2) ) D x(2) tto(t (2) ,x (2) = 0) 

J — oo J -co 

+ 7 2 / dt' / dt^ / dt^ / dt^ 

J —oo J —oo J —oo J— oo 

xD^GMk^x;?^' = 0) G osc (t';t {2) ) D^D^G^t^, x {2) = 0;t (3) ,x (3) = 0) 
xG osc (t (3) ;t (4) )D xW u (t^,x^ =0) 
+ ... (152) 



dt 2 



+ uo 2 {t) 



u(t, x = 0) = 7 



du 
dx 



d 2 



Ui(t, x = 0) = 7 



<9x 



d 2 



u n+1 (t,x = 0) = 7 



9a; 
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This expansion contains the infinite sum 

G sc G osc 

+ lG 0SC {D x iD x nGbuik)G sc 

+ l 2 G osc (D x jD x jiGbuik)G osc (D x jD xI iGbuik)G osc 

+ l Z G osc (D xI D xII Gbuik)G osc (D xI D xII Gbuik)G osc (D xI D xII Gbuik)G osc + . . . (153) 

Here the notation D x j and D x u indicates that the spatial derivative acts on the spatial 
component of the first and the second argument of G^uik, respectively. We may decompose 
the internal bulk propagator lines into a singular and a regular parts 

D^uG^^K^ + K^. (154) 

Consider for example the massless propagator 

[dt -d 2 x ]-\t,x-t',x') = (1/2) 6(t-t r ) 6((t-t') 2 -(x-x') 2 ), (155) 

or the uniform mass propagator 

[d 2 -d 2 x + m 2 ]-\t,x-t',x') = 

(1/2) 6{t - t') e(y/{t - Vf -{x- x') 2 )j Q (m^(t - t') 2 - (x - x'f) . (156) 

These are the propagators for the entire infinite plane 1Z 2 . We form the retarded Dirichlet 
propagator by using the method of images so that 

G D (t, x- 1\ x') = Goo(t, x- 1', x') - Goo(t, x; t', -x'). (157) 

It turns out that the singular (local) part has the form 

K { S\t-t') = -8'{t-t'), (158) 

no matter what the propagator is, that is, irrespective of the particular form of the function 
m 2 (x) in eqn. (|141j) . By contrast, depends sensitively on the particular form of m 2 (x) 

and encodes all the non-locality arising from propagation from the brane into the bulk and 
back again. 
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+ 5 + L\ + ... 

i r 





FIG. 4: Perturbation expansion for the full brane-brane two-point function as a power 
series in the brane-bulk coupling 7. The first line shows three increasingly complete forms 
for the brane-brane propagator. The most complete is the full propagator, denoted by the triple 
line, including both dissipative and nonlocal effects due to the interaction with the bulk. The 
double line includes only the dissipative effects of the bulk interaction, and the single line does not 
include any correction at all from the brane-bulk interactions. For the simplest case, the single 
line represents (D% + an d the double line (Df + jD t + Wq)^ 1 . More generally, ujq and 7 may 

depend on time. On the second line we show the diagrammatic decomposition of the full bulk-bulk 
propagator K (with both endpoints on the brane and the brane-bulk coupling factors included), 
represented by the dashed curve. This propagator decomposes into its purely local singular part, 
represented by the cross, and its nonlocal, regular part, represented by the dotted line. The third 
line illustrates the perturbative expansion for the full brane propagator. On the fourth line the 
partial sum over the local insertions is carried out. The last line shows the expansion of the full 
propagator in terms of the dissipation corrected propagator and the nonlocal corrections. 
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We use the massless propagator 



Gd(x,x') = GWx,x') - Goofex'^) 



(159) 



because the correction due to mass or, more generally, any potential contributes only to 
the nonlocal part of Kbuik where t and t' do not coincide. 
We have 



K(t,t') = d x d x ,G D {t,x-t',x' 



— lim 

2 a;^0+ 



x=x'=0 



= i \im + lim_ d x d x , (o(t - t') ■ l^6(t - t' + x - x') - 9(t - t' - x - x') j^j 



limj -5'{t -t' + x-x')- 5'{t - t' - x - x') 



= -5'{t-t'). 



(160) 



The order of the limits is arbitrary; however, by choosing a certain order as we have done 
here, we reduce the number of terms to consider. 

We now consider the contribution for t > t' where we have used the massive propagator 



where 



G 0O (x,^) = -9(t-lf) 0(r)Jo(mr), 



r = y/(t-t') 2 -(z-z') 2 . 



(161) 



(162) 



One has 



\_d__d_ 

2 dx dx' 



= m 



Jg^m(t — t 
mJ' (m(t — t' 



J (m^/ (t - t') 2 - (x - x') 2 ) - J (m^ (t - 1') 2 - (x + x') 2 ) 
f<\ 9 9 



x=x'=0 



(163) 



(t - t>) 

We observe that the regular part here is in fact regular as t — > t'+, since J' Q {z) ~ z. This is 
a universal feature. 
Since 



G, 



D 2 + co 2 (ty 



(164) 
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it follows that 

G osc = G osc + G osc (^K^ 9 )G osc + Gosc^K^^jGosc^yK^. 9 JG osc + . . . 

= A 2 + 7A + ^) (165) 

and 

G osc = G osc + G^K^G^ + G^K^G^K^G^ + .... (166) 

The above expansions may be interpreted pictorially by means of Feynman diagrams, 
as indicated in Fig. HJ The propagator for the oscillator with the coupling to the bulk 
ignored given in eqn. (|164|) is denoted by a single vertical line localized on the brane. The 
dashed curves represent propagation of bulk gravitons and their coupling to the brane. We 
decompose the bulk propagator into two parts: a purely local singular part, denoted by 
a cross insertion, and a regular, nonlocal part denoted by a dotted curve, as indicated on 
the second line of the figure. The third line shows how the full two-point function on the 
brane (indicated by the triple line) may be expressed as an infinite sum. The fourth line 
indicates a partial infinite sum with only the local insertions included yielding an improved 
propagator on the brane (denoted by by the double vertical line). In terms of the effective 
equation on the brane, this infinite summation simply represents the addition of a local 
dissipation term to the oscillator equation of motion. Finally, on the last line, an infinite 
sum adding the nonlocality is taken to calculate the exact full two-point function on the 
brane. 



V. GENERALIZATION TO MOVING BOUNDARIES 



For computing braneworld cosmological perturbations, we are interested in systems with 
a moving boundary. In the picture where Birkhoff 's theorem is exploited so that the bulk is 
stationary and the brane moves, except for the boundary condition imposed by the moving 
brane, the propagation in the bulk is trivial. One simply propagates the bulk gravitons 
using the coordinate system that best exploits the symmetries of AdS. 

For the scalar wave equations considered here, the closest equation to the case of AdS 
gravitons is 

dt 2 dx 2 x 1 



u(t, x) = 0, 



(167) 
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where k corresponds to the magnitude of the wave number in the three transverse spa- 
tial dimensions. The boundary is described by a trajectory (tb), which may be 
parameterized by proper time r where 

^ = a(x b ) • y/l-x', 2 (168) 

and a(xb) is the conformal scalar factor of the metric. We also define a barred proper time 
f where this conformal factor has been omitted. The oscillator is treated as in the previous 
section except that now the proper time r takes the place of ordinary time. We parameterize 
x(t) = (tb(r), Xb(r)j . Most everything discussed in the previous section for the case of 
a stationary boundary generalizes straightforwardly to the case of a moving boundary. 
The only slight complication is the computation of the Dirichlet Green's functions and 
the behavior of the Ui n and u out solutions, which are represented by the external lines in 
our diagrammatic expansion. However, this may be accomplished by placing a virtual 
source an infinitesimal distance just beyond the boundary. This method, which is a sort of 
generalization of the method of images, is originally due to d'Alembert.|lll] 

For the external lines in our diagrammatic expansion in the oscillator/bulk coupling 
strength 7, we seek two complete bases of "in" and "out" states, characterized by the 
asymptotic behavior on past and future null infinity, respectively. The solutions satisfy 

d 2 d 2 



CbulkU(t,x) 



u(t,x) = (169) 



dt 2 dx 2 

where u can be either Ui n or u ou t and satisfy the hard boundary condition 

C-BMrdV = (170) 

on the boundary x = x_b(f ). 

By hard we mean the condition that is local on the boundary that agrees with the exact 
boundary condition (which we shall call "soft") in the short-distance or high frequency 
limit. For the examples considered in this paper this is a Dirichlet boundary condition, 
where u = on the boundary, but Neumann and mixed boundary conditions are allowed, 
because these do not relate the behavior at different points on the boundary in a nonlocal 
way. 

For a stationary boundary the asymptotic states u ou t and Ui n may be readily constructed 
by separation of variables. For a moving boundary accounting for multiple reflections is 
less trivial. 
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FIG. 5: Method of virtual sources. We illustrate how the bulk propagator from x' to x may be 
calculated for a moving boundary using the method of virtual sources. In general we do not know 
the form of the propagator satisfying the relevant boundary conditions. Therefore, to compute 
G(x, x')> we start with an arbitrary propagator, in general violating the boundary conditions, and 
then place a virtual source just infinitesimally beyond the boundary to correct for this violation. 
Here a corrective virtual source is required on that portion of the boundary intersecting with 
both the forward lightcone of x' and the past lightcone of x, which we parameterize by f. This 
virtual source (for the case of Dirichlet boundary conditions) is a dipole layer, whose strength as 
a function of position on the boundary is determined by solving a Volterra integral equation of 
the second kind. 

Before calculating the extensions of u in and u out based on their behaviors in the neigh- 
borhood and past and future null infinity, respectively, we first construct the retarded 
Green's function satisfying the Dirichlet boundary condition on the null boundary 



starting from the retarded Green's function that would hold if there were no boundary 



We proceed by postulating a virtual source just beyond the boundary, propagated ac- 
cording to G^, which serves to enforce the boundary conditions on the moving boundary. 
For Dirichlet boundary conditions this is a dipole layer of strength characterized by the 
linear density m(f) where the dipole is inward and normally directed. The field on the 
boundary of a pointlike dipole of unit strength at r' just beyond the boundary is of the 




(171) 



(172) 
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form 

U{f-f')+6{f-f')R{f,f') (173) 

where R(t, t') is a smooth, regular function without any divergence at f = f'. [For the case 
of Neumann or mixed boundary conditions, we would use a monopole rather than a dipole 
layer.] 

We write as an ansatz 

G D (x,x') = Goofex') + C df (n(f) • V 2 ) G 00 (x,x(r)) m(f;x,x') (174) 

as indicated in Fig. Here n(r) is the inward normal vector on the boundary and V2 
denotes derivative with respect to the second argument of the Green's function. Here the 
interval [r a , TjJ, which can be empty for certain choices of x, x', is the intersection of the 
boundary worldline, the interior of the forward lightcone of x', and the interior of the past 
lightcone of x. The virtual source m(f; x, x') is chosen so that the boundary conditions are 
satisfied. This is accomplished by solving the integral equation 

^m(r;x,x') = Goo (x(f) , x') + [ df' (n(f) ■ V 2 ) Gj^f), x(t')) m(f';x,x'). (175) 

It is understood that the singular part, which occurs on the left-hand side of the equation, 
has been removed from the integral kernel. This is a Volterra inte gral equation of the 
second kind, which is always guaranteed to have a unique solution 12j. The "in" and 
"out" functions may be constructed analogously using an integral equation. 



VI. DISCUSSION 



One of the most interesting questions of braneworld cosmology is the nature of the 
cosmological perturbations predicted and how these might differ from those possible in 
models not having a large extra dimension. In order to carry out the required calculation, 
the brane-bulk interaction must be taken into account, preferably without resorting to ap- 
proximations based on lower- dimensional effective descriptions. The techniques developed 
in this paper, illustrated using scalar toy models in which a (1+1) -dimensional scalar field 
theory is coupled to a finite number of degrees of freedom on its boundary, constitute a 
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step toward achieving such a complete calculation and demonstrate some of the qualitative 
effects that one might expect. 

In section II we studied some simple time-independent systems. It was seen how coupling 
the boundary (brane) degrees of freedom to the continuum (bulk) gave rise to dissipative 
effects from a purely brane point of view. In addition, lack of spatial homogeneity in the 
continuum (as would result from the "warp factor" of Randall- Sundrum- like cosmologies) 
led to nonlocal effects, which can be encoded into a spectral density p{uS). In section III we 
saw how new effects could arise in time-dependent systems (as would arise in an expanding, 
rather than static Minkowski, geometry on the brane). The discrete degrees of freedom 
on the brane, rather than being coupled at the same strength infinitely far into the past, 
could initially be uncoupled, interacting effectively for only a finite amount of time, so that 
a sort of S matrix linking the "in" oscillator states to the "out" oscillator states for the 
brane and bulk dynamics would connect observables today with initial conditions for the 
brane and the bulk. When we observe the braneworld perturbations today, for the simplest 
case where these are Gaussian, it suffices to calculate the expectation values of operators 
quadratic in a branefiUt and a braneout . [Note that our simplified notation here assumes that 
the wave number in the three transverse spatial dimensions k has been fixed and ignores 
the fact that these operators also have an integer index to account for the possibility of 
several modes on the brane. Similarly, the bulk operators have a discrete index to account 
for the several bulk graviton polarizations and other possible modes propagating in the 
bulk. Since the generalization is straightforward, we suppress these complication in the 
notation that follows.] In terms of the "in" state, a branefiUt and a\ rane out may be expressed 
as a linear combination of "in" operators, A brane>in and A bu i k>ini normalized so that 

[A brane , in , Al ranein ] = 1, (176) 

and 

[Abuik,in,Al ulkin ] = 1, (177) 

where A brane)in is constructed entirely as a linear combination of a brane ^ n and a brane in , and 
likewise A bu ik t i n is constructed entirely as a linear combination of a bu ik,i n (u) and al ulkiin (u), 
where u > 0. There are three possibilities: firstly, 

&brane,out = COS 9 A brane j n + sin 9 A bu ikj n (178) 
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where < 9 < n/2; secondly, 

a br ane,out = cosh f A braneM + sinh £ Al ulkM (179) 
where < £ < +oo; and thirdly, 

&brane,out 

sinh f Aj mne in + cosh f (180) 

where < £ < +oo. Physically, when we measure the perturbations today, we observe in 
part the character of the initial state on the brane and in part that of the initial state of 
the bulk. Their relative importance can be read off from this S-matrix. 

In section IV we dealt with reflections in the bulk for the time-dependent situation. Here 
a wave initially propagating away from the brane is scattered, reflected, or diffracted— 
depending on which term one prefers — by the bulk so as to propagate back toward the 
brane, and vice versa. A perturbation expansion was developed where "hard" boundary 
conditions are imposed at zeroth order and the "soft" corrections are included by successive 
approximations. A "hard" boundary condition is defined as one that is completely local 
on the boundary, such as the Neumann or Dirichlet boundary condition, or a mixture 
of the two. These are "non-dynamical" because they do not interrelate the boundary 
behavior at distinct points on the boundary. By contrast, "soft" boundary conditions are 
dynamical and nonlocal, resulting when the degrees of freedom on the boundary cede to 
the incoming waves. Asymptotically, away from the resonances of the degrees of freedom 
on the boundary, the exact "soft" boundary conditions approach a corresponding "hard" 
boundary condition, which we use at zeroth order in our perturbation expansion. 

In section V we indicated how to deal with the complication of a moving boundary. In 
most respects, the techniques of the previous sections generalize straightforwardly. The 
time for the degrees of freedom on the boundary becomes proper time. The principal 
difficulty in carrying out the generalization to a moving boundary having an arbitrary 
trajectory (corresponding to an arbitrary expansion history for braneworld cosmology) is 
the calculation of the bulk Green's functions satisfying the hard boundary conditions. This 
is accomplished by the method of virtual sources, where an arbitrary Green's function (in 
general not respecting the boundary condition) is used and then a virtual source placed 
just infinitesimally behind the boundary is postulated in order to emit into the bulk waves 
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that correct for the violation of the boundary condition due to using the wrong Green's 
function. The source that achieves such a correction is calculated by solving an integral 
equation of the Volterra form of the second kind. 

We are in the process of generalizing this work to include gravitational gauge fixing so 
that it can be applied to calculating braneworld cosmological perturbations to linear order. 
This work will be the subject of a future publication. 
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